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ABSTRACT 

We investigate the effects of the Cosmic Microwave Background (CMB) radiation 
held on the collapse of prestellar clouds. Using a semi-analytic model to follow the 
thermal evolution of clouds with varying initial metallicities and dust contents at 
different redshifts, we study self-consistently the response of the mean Jeans mass at 
cloud fragmentation to metal line-cooling, dust-cooling and the CMB. 

In the absence of dust grains, at redshifts z ^ 10 moderate characteristic masses 
(of 10s of M Q ) are formed when the metallicity is 1O _4 Z Z 10" 2 5 Z Q ; at 
higher metallicities, the CMB inhibits fragmentation and only very large masses (of 
~ 100s of Mq) are formed. These effects become even more dramatic at z > 10 and 
the fragmentation mass scales are always ^ 100s of Mq, independent of the initial 
metallicity. 

When dust grains are present, sub-solar mass fragments are formed at any redshift 
for metallicities Z ^ 10~ 6 Z Q because dust-cooling remains relatively insensitive to the 
presence of the CMB. When Z > 10~ 3 Z Q , heating of dust grains by the CMB at z ^ 5 
favors the formation of larger masses, which become super-solar when Z ^ 10~ 2 Z Q 
and z ^ 10. Finally, we discuss the implications of our result for the interpretation of 
the observed abundance patterns of very metal-poor stars in the galactic halo. 

Key words: Cosmology: cosmic microwave background, galaxies: evolution, stars: 
formation, Population II, ISM: abundances, dust 



1 INTRODUCTION sets the mean Jeans mass at cloud fragmentation (see Bon- 

nell, Larson & Zinnecker 2007 for a comprehensive review), 

The origin of the stellar Initial Mass Function (IMF) and its 

variation with cosmic time or with diverse environmental / Tf r ,„ \ 3 / 2 / rtr raa \~ 1 / 2 

M frag = ^~7OOM (^!) (lo^) (D 



conditions still lack a complete physical interpretation. Ob- 
servationally, the present-day stellar IMF appears to have 
an almost universal profile, characterized by a power-law 



where the Jeans length is defined as, 



at large masses and flattening below a characteristic mass / -xkT ^ ^ 2 

of Af c h ~ IMq with a plateau which extends down to ^ J ~ I ^ mH Qp 
the brown dwarf regime (see Chabrier 2003 and references y 
therein). This mass scale, which represents the "typical" The process of gravitational fragmentation can be investi- 
outcome of the star formation process, appears to be re- gated using thermal physics within the framework of the 
markably uniform in diverse environments where the ini- classical Jeans criterion: fragmentation occurs at the end 
tial conditions for star formation can vary considerably (see of an efficient cooling phase, when the temperature de- 
Elmegreen, Klessen & Wilson 2008). creases for increasing density and the effective adiabatic in- 

Among the many proposed explanations, the origin of dex switches from 7 < 1 to 7 > 1 (the pressure p a p 7 , 
the characteristic stellar mass and the broad peak of the IMF where p is the density) ; the typical fragmentation mass-scale 
are best attributed to gravitational fragmentation, which is given by the thermal Jeans mass at that epoch (Larson 

2005; Schneider et al. 2002). The appearance of a charac- 
teristic mass at the inflection point in the equation of state 
has been also demonstrated numerically in 3D simulations 
* E-mail: raffa@arcetri.astro.it by Jappsen et al. (2005), supporting the idea that the dis- 
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tribution of stellar masses depends mainly on the thermo- 
dynamical state of the star-forming gas. 

In the physical conditions which apply to star forming 
regions in the present-day Universe, the application of the 
above criterion leads to a characteristic fragmentation mass 
of Mfrag ~ 2Mq. About a half of this gas is expected to 
be ejected by a protostellar outflow without accreting (e.g. 
Matzner & McKee 2000). 

If we apply the same criterion to primordial star form- 
ing regions we find a fragmentation mass of order Mf rag ~ 
7OOM , which corresponds to the Jeans mass at the tem- 
perature and density where cooling by molecular hydrogen 
becomes less efficient due to the NLTE-LTE transition in 
the level populations of the molecule. These clumps are the 
immediate progenitors of the stars, in agreement with the 
results found by 3D numerical simulations (Abel, Bryan 
& Norman 2002; Bromm, Coppi & Larson 2002; Yoshida, 
Omukai & Hernquist 2008). 

Thus, between the primordial and the present-day con- 
ditions there is a clear transition in fragmentation scales 
which is believed to be mainly driven by the progressive en- 
richment of star forming clouds with metals (Bromm et al. 
1999; Omukai 2000; Schneider et al. 2002; Bromm & Loeb 
2003; Santoro & Shull 2006; Smith, Sigurdsson & Abel 2008) 
and dust grains (Schneider et al. 2003; Omukai et al. 2005; 
Tsuribe & Omukai 2006; Clark, Glover & Klessen 2008). 

Note, however, that the presence of gas-phase metals 
can not by itself lead to solar or sub-solar characteristic 
masses. In fact, metal line-cooling is efficient in the ini- 
tial evolution of collapsing gas clouds, for densities in the 
range 10 4 - 10 6 cm , and induces fragmentation into rela- 
tively large clumps, with characteristic fragment masses of 
a few 10s of Mq up to initial gas metallicities Z <, 10~ 2 Zq 

(Schneider et al. 2006). Conversely, the formation of solar or 
sub-solar mass fragments might operate already at metallici- 
ties Z CT = 10~ 6 Zq in the presence of dust grains, which pro- 
vide an efficient source of cooling at densities £ 10 13 cm~ 3 . 

Since dust can be promptly synthesized in the ejecta of 
primordial supernovae (Todini & Ferrara 2001; Nozawa et 
al. 2003; Schneider, Ferrara & Salvaterra 2004; Bianchi & 
Schneider 2007), it is likely that its contribution to ther- 
mal and fragmentation history of collapsing clouds at high 
redshifts had been relevant. 

An additional redshift-dependent effect is induced by 
the Cosmic Microwave Background (CMB) which sets an ef- 
fective temperature floor increasing with redshift. Clarke & 
Bromm (2003) have proposed a simple formula for the char- 
acteristic stellar mass which is identified by the Jeans mass 
fixed jointly by the gas temperature and the pressure set by 
the weight of the overlying baryons as the gas collapses in 
the parent dark matter halo. In this simple model, they do 
not follow the gas thermal evolution and the temperature is 
taken to be the maximum value between the CMB floor and 
either 200 K for Hb-cooling or 10 K for CO-cooling. In the 
latter case, the characteristic stellar mass is set by the CMB 
temperature at all redshifts 2 > 2.7 and it is predicted to be 
> 10M Q at z > 10. 

The CMB influences the thermal evolution of collapsing 
gas clouds in two possible ways: (i) it affects the atomic and 
molecular line level populations (ii) it heats the dust grains. 
In Omukai et al. (2005), we have estimated the impact of the 



CMB temperature floor at z = 20 on the evolution of prestel- 
lar clouds with varying initial metallicity. It was shown that 
the CMB affects the cloud evolution at low temperatures, 
which are reached for relatively large values of the initial 
metallicity, Z > 1O~ 3 Z . 

Smith et al. (2009) have investigated the effects of the 
CMB on the collapse of prestellar clouds with varying initial 
metallicities by means of numerical simulations. They con- 
firm that due to metal line-cooling, the gas becomes ther- 
mally coupled to the CMB for Z > 10~ 2 !i Z e early in the 
evolution of the clouds when the density is still low. Note 
that a similar effect has been recently found by Jappsen 
et al. (2009) using 3D simulations starting from different 
initial conditions, namely protogalaxies forming within a 
previously ionized HII region that has not yet had time to 
cool and recombine. Even in the latter models, the gas cools 
rapidly to the CMB floor and no fragmentation occurs dur- 
ing the collapse for metallicities as high as O.1Z . 

Smith et al. (2009) conclude that the CMB inhibits frag- 
mentation at large metallicities and that there were three 
distinct modes of star formation at high redshift: a primor- 
dial mode, producing very massive stars at very low metal- 
licities [Z < 10 -3 ' 75 Zq); a CMB-regulated mode, produc- 
ing moderate mass (10s of Mq) stars at high metallicites 
(Z > 10 -2 ' 5 Zq at redshift 2 ~ 15 — 20); and a low-mass (a 
few Mq) mode existing between those two metallicities. As 
the Universe ages and the CMB temperature decreases, the 
range of the low-mass mode extends to higher metallicities, 
eventually becoming the only mode of star formation. 

In their simulations, however, Smith et al. (2009) con- 
sider the thermal evolution in the presence of molecular 
coolants (H2 and HD) and gas-phase metals only. They do 
not include H2 formation heating, which affects the thermal 
evolution at moderate densities > 10 s cm~ 3 in the primordial 
gas and at lower densities in the metal-enriched gas. More- 
over, they do not take into consideration the possible pres- 
ence of dust grains which significantly affect the gas thermal 
evolution even at very low metallicities via cooling by ther- 
mal emission and efficient H2 formation on their surfaces 
(Schneider et al. 2003; Omukai et al. 2005; Schneider et al. 
2006). 

From an observational point of view, the number and 
properties of very metal-poor (VMP) stars in the Galactic 
halo with [Fe/H] < —2 seem to require a higher charac- 
teristic stellar mass in the past, with a redshift-modulation 
consistent with that expected by the CMB (Hernandez & 
Ferrara 2001). Recent analyses of the surface elemental com- 
position of very metal-poor halo stars (Cayrel et al. 2004) 
show that the observed abundance pattern appears consis- 
tent with the predicted yields of Pop III stars with char- 
acteristic mass of « IOMq (Heger & Woosley 2008). Sim- 
ilar conclusions have been derived analysing the fraction 
of carbon-enhanced stars (CEMP, Tumlinson 2007; Komiya 
et al 2007). CEMP stars are a subset of metal-poor stars 
that show enhanced carbon-to-iron abundances ([C/Fe] > 1) 
(Beers & Christlieb 2005). These stars are thought to be 
metal-poor low-mass stars (< O.8M ) that have acquired C 
enhancements at their convective surfaces by capturing the 
C-rich ejecta of an AGB companion (1.5M < M* < 8M ). 
Because the binary system that produces a CEMP star re- 
quires both a low-mass star and an intermediate mass star, 
the incidence and chemical abundance signatures of CEMP 
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stars reflects the underlying IMF in the range 1M© < M* < 
8M Q (Lucatello et al. 2005; Tumlinson 2007; Komiya et al. 
2007). 

On the basis of these considerations, Tumlinson (2007) 
suggests that the CMB is shaping the characteristic mass of 
the IMF and proposes the following parametrization, 



M ch /M e = 1.06M Q 



max[(l + z)2.73 K,8 K] 
10K 



(2) 



which implies that M ch = 0.72, 6.87 and 20.63M© at z = 
0, 10 and 20. 

The question is whether this observationally motivated 
CMB-regulated IMF can be reconciled with theoretical mod- 
els for the thermal evolution and fragmentation of star- 
forming clouds at very low metallicities. 

Here we intend to revisit the role of the CMB in the evo- 
lution of protostellar clouds investigating in a self-consistent 
way the response of the fragmentation mass scale to metal 
line-cooling, dust cooling, and the CMB. 

The paper is organized as follows: in section [2] we 
present the model adopted to follow the thermal and chem- 
ical evolution of collapsing protostellar clouds, in section [3] 
we illustrate the main results of the analysis, and in section 
[4] we discuss their implications. 



2 THE MODEL 

The thermal and chemical evolution of star forming gas 
clouds is studied using the model developed by Omukai 
(2000) and improved in Omukai et al. (2005) . The collapsing 
gas clouds are described by a one-zone approach where all 
physical quantities are evaluated at the center as a function 
of the central density of hydrogen nuclei, The tempera- 
ture evolution is computed by solving the energy equation, 

de d 1 , . 

d-t = ~ P di-p- A ^ (3) 

where the pressure, p, and the specific thermal energy, e, are 
given by, 



P = 



pkT 



1 



kT 



7ad - ljU 7BH ' 



(4) 



(5) 



and p is the central density, T is the temperature, 7 a d is the 
adiabatic exponent, p is the mean molecular weight and ran 
is the mass of hydrogen nuclei. The terms on the right-hand 
side of the energy equation are the compressional heating 
rate, 



dp 
~dt 



P_ 

tft 



with 



ts 



3tt 
32Gp' 



(6) 



and the net cooling rate, A nct , which consists of three com- 
ponents, 



An ne + Ao 



(7) 



The first component, Aii ne , represents the cooling rate due 
to the emission of line radiation, which includes molecular 



line emission of Ha, HD, OH, H2O and CO, and atomic fine- 
structure line emission of CI, CII and OI. Following Omukai 
et al. (2005), H2 collisional transition rates and HD parame- 
ters are taken from Galli & Palla (1998). The second compo- 
nent, A cont , represents the cooling rate due to the emission 
of continuum radiation, which includes thermal emission by 
dust grains and H2 collision-induced emission. The last term, 
A c hem, indicates the cooling/heating rate due to chemical 
reactions. When the gas cloud is optically thick to a spe- 
cific cooling radiation, the cooling rate is correspondingly 
reduced by multiplying by the photon escape probability. 
Unless otherwise stated, the treatment of these processes 
are the same as in Omukai et al. (2005), to which we refer 
the interested reader for further details. 

In particular, we adopt the reduced chemical model for 
low-metallicity clouds described in Section 3.2 of Omukai et 
al. (2005). This model reproduces the temperature evolu- 
tion and the abundances of the dominant chemical species 
following a network of 54 reactions (instead of the full 500 
reactions of the original model) for 26 species (instead of 
50). 

In the present analysis, we are interested in investi- 
gating the effects of the CMB in the thermal evolution of 
the clouds at different metallicities. The CMB enters in our 
model in two ways: first through the level population of 
atomic /molecular lines, and second through the heating of 
dust grains. These effects prevent the gas and dust temper- 
ature from falling below the radiation temperature. In addi- 
tion, the higher dust temperature reduces the H2 formation 
rate by surface reactions. The former effect is included in 
calculating the level populations of atomic and molecular 
species as in Omukai (2001, see Appendix B). The latter 
effects are included through the energy-balance equation of 
grains, 



4 (J Tcry ^SIT -A-ff 



^dust + 4 a T ra( j Av gr 



(8) 



where the left-hand side is the cooling rate of grains and 
the right-hand side is the heating rate by collisions with gas 
particles and by the CMB (see Omukai et al. 2005), and we 
have made the approximation that the Planck mean opacity 
for absorption is the same as for emission. 

For gas at solar metallicity, the abundances of carbon 
and oxygen are yc = 0.927 x 10~ 4 and yo = 3.568 x 10" 4 . 
The mass fraction of dust grains is 0.939 x 10~ 2 below 
the water-ice evaporation temperature, and it decreases for 
higher temperatures as each grain component evaporates 
(Pollack et al. 1994). For lower metallicity gas, these val- 
ues are reduced proportionally and we denote the corre- 
sponding abundances relative to the solar value as, [M/H] = 
log 10 (Z/Zo). Strictly speaking, we do not expect the re- 
gions enriched by the first supernovae to have a solar abun- 
dance pattern and the same dust grains present in the lo- 
cal neighborhood. Yet, a more detailed analysis which self- 
consistently accounts for the elemental abundances and dust 
grain properties produced by metal-free supernovae with 
masses 12 - 30 M Q and 140 - 260 M shows that the ther- 
mal evolution of protostellar coulds is qualitatively the same 
(see Schneider et al. 2006). 
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Figure 1. Thermal evolution as a function of the central gas number density for protostellar coulds with initial metallicities [M/H] = 
-5, -4.5, -3, -2, -1 (from the top-left to the bottom-right). Temperature divided by mean molecular weight is shown. For our assumed 
helium abundance (j/He = 0.0833), the mean molecular weight changes from fj, = 1.23 (fully atomic hydrogen) to 2.29 (fully H2). In each 
panel, the different curves correspond to the evolution of clouds at 2 = (solid), 5 (dotted), 10 (short-dashed), 20 (long-dashed), and 30 
(dot-dashed). The points identify the epochs of fragmentation induced by line-cooling (dots) and dust-cooling (triangles). The diagonal 
dotted lines identify values of the Jeans masses corresponding to given thermal states. 
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3 RESULTS 

In this section we discuss the model predictions for the ther- 
mal evolution of low-metallicity protostellar clouds at dif- 
ferent redshifts, i.e. exposed to a CMB radiation field at 
different temperatures. We then discuss the resulting char- 
acteristic fragmentation mass scales. 

3.1 Thermal evolution 

In Fig. [l]we show the thermal evolution of collapsing proto- 
stellar clouds as a function of the central gas number density. 
The temperature divided by the mean molecular weight T/fi 
is shown as its variation with density is directly related to 
the effective adiabatic index 7 — 1 = dln(T/^)/dlnnH. Each 
panel represents the evolution for a specific initial metal- 
licity, ranging from [M/H] = -5 (top-left) to [M/H] = -1 
(bottom-right). For each metallicity, we have explored the 
variation with redshift and the different curves in each panel 
represent clouds collapsing at z = (solid), 5 (dotted), 
10 (short-dashed), 20 (long-dashed), and 30 (dot-dashed). 
Thus, the effects of the CMB can be quantified by the devi- 
ations of the curves from the reference solid curve because 
at z — the CMB temperature is too low to have any ap- 
preciable consequence on the clouds thermal evolution. 

As expected, the effects of the CMB increase with red- 
shift (i.e. higher radiation temperature) and metallicity (i.e. 
lower gas temperature). However, even at metallicities as 
small as [M/H] = —5, the thermal evolution at z = 30 shows 
deviations from the z — track, despite the fact that at this 
low metallicity the gas temperature is higher than the CMB 
value. The deviation is caused by the suppression of H2 for- 
mation on dust grains due to the increased dust tempera- 
ture, which is almost equal to the CMB temperature. The 
H2 formation rate coefficient on dust grains is expressed as 
(Tielens & Hollenbach 1985), 

_ 6 x lO- 17 cm 3 s- 1 (r/3OOK) 1/2 (Z/Z )/ a 

H2,gr ~~ 1 + 4 x 10" 2 (T + T gr ) 1/2 + 2 x 10" 3 T + 8 x 10" 6 T 2 ' 

with 



Ja 1 +exp[7.5 x 10 2 (l/75 - 1/T gr )]' 

Thus, when T gr > 75 K, / a becomes very small. If the H2 
formation rate is not reduced at high dust temperatures, as 
suggested by Cazaux & Tielens (2004), the dependence of 
the thermal evolution on the CMB at these low metallicities 
would disappear. 

For initial metallicities [M/H] < —3, the effects of the 
CMB on the cloud thermal evolution are limited to densities 
71h < 10 8 cm" 3 , in the regime where line-cooling dominates. 
At higher metallicities, the minimum temperature induced 
by dust cooling reaches the CMB value. Thus, the deviations 
from the z = track are more pronounced and are present 
at all densities riH 10 3 cm -3 , up to the regime where cool- 
ing is dominated by dust grains. Note that at [M/H] ^ —2 
and z > 20, decrease in T/^i beyond the point where the 
temperature reaches the CMB floor is due to the increase in 
)i as a result of H2 formation. 

In each panel, the points identify the epochs of fragmen- 
tation induced by line-cooling (dots) and by dust-cooling 
(triangles), while the diagonal dotted lines represent some 




number density log n H (cm 3 ) 

Figure 2. Evolution of the adiabatic index 7 for clouds with 
initial metallicities [M/H] = —5, —4, —3, —2, —1 (from top to bot- 
tom). In each panel, the different curves show the evolution for a 
fixed metallicity and redshift z = (solid), 5 (dotted), 10 (short- 
dashed), 20 (long-dashed), and 30 (dot-dashed). The horizontal 
dotted lines indicate the values limiting the fragmentation epoch, 
i.e. a cooling phase with 7 < 0.8 followed by the condition that 
7 > 0.97 (see text). 

indicative values of the Jeans mass corresponding to the 
thermal states. As expected on the basis of arguments given 
in section 1, the fragmentation epochs correspond to the 
inflection points of the equation of state. 

More quantitatively, we identify the fragmentation 
epochs by the condition that the adiabatic index becomes 
7 > 0.97 after a phase of cooling (where 7 << 1). The choice 
of a threshold value slightly smaller than unity is motivated 
by the fact that when the gas reaches the CMB tempera- 
ture, 7 may approach unity asymptotically, but never ex- 
ceeds it. As an additional condition for fragmentation, we 
impose that 7 < 0.8 during the cooling phase preceding the 
fragmentation epoch: this is to eliminate "false" fragmenta- 
tion points where 7 is less than 1 only for a short period of 
time. The evolution of the adiabatic index 7 as a function of 
the central number density for the same set of models shown 
in Fig. [T] is shown in Fig. [2] 

For all but the lowest (highest) metallicity models, 
the two fragmentation phases (line-induced, at low den- 
sities, and dust-induced at high-densities) are preserved 
even in the presence of the CMB. When [M/H] = —5 and 
z — 30, however, there appears to be an intermediate frag- 
mentation phase which is due to H2O cooling at densities 
10 7 - 10 8 cm . At lower redshifts/higher metallicities, this 
phase is no longer present because at these intermediate den- 
sitites the temperature is already too low (due to H2 /metal 
line-cooling) for H2O cooling to have any appreciable ef- 
fect. However, as it can be inferred from the behaviour of 
the dot-dashed line in the top panel of Fig. [2] the interme- 
diate fragmentation phase is most probably artificial as it 
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depends very much on the fragmentation condition. With a 
stricter condition on the duration of the cooling phase, such 
as 7 < 0.7, it would not be classified as a fragmentation 
phase. 

Finally, as it can be inferred from the bottom panels of 
Figs[T]and[2] when [M/H] > -2 and z > 10, the line-induced 
and dust-induced fragmentation phases merge and the ther- 
mal evolution is controlled by the CMB temperature. 

3.2 Fragmentation mass 

In the present context, where we assume that thermal pres- 
sure is the main force opposing gravity, the characteristic 
mass of each fragment is given by the Jeans mass at fragmen- 
tation epoch, given in eq.[T] Fig. [3] shows the fragmentation 
masses as a function of metallicity and redshift. In the left 
panel, the upper and lower branches are associated to line- 
induced and dust-induced fragmentation, respectively. Since 
line-induced fragmentation occurs at relatively low densities, 
the corresponding fragment masses are large. In particular, 
at the lowest metallicities we find that, at any redshift, the 
characteristic fragmentation mass scale, induced by H2 line 
cooling, is Mf rag ~ 10 3 Mq, consistent with the results of 
numerical models. Still, even for metallicities in the range 
—5 ^ [M/H] ^ —2, metal-line cooling leads to fragment 
masses Mf rag ^ 50M© (Schneider et al. 2006). Dust-induced 
fragmentation, which occurs at much higher densities, leads 
to significantly lower mass scales. 

The CMB affects the amplitude of the characteristic 
mass in both line- and dust-induced phases. At any metal- 
licity in the range —5 [M/H] ^ — 1, the presence of 
the CMB leads to larger characteristic masses in the line- 
induced phase, with a deviation from the reference 2 — 
track which increases with metallicity and redshift. Even 
for clouds with metallicity [M/H] = -4, the characteristic 
mass is more than a factor 2 larger at 2 ^ 20. As it has 
already been discussed, the main effect of the CMB at these 
low metallicities is the suppression of H2 formation on dust 
grains. In the dust-induced phase, instead, the effects of the 
CMB on the characteristic mass start to be appreciable only 
when the metallicity is [M/H] —3, when the temperature 
dip due to dust cooling hits the CMB floor. At [M/H] = -3, 
the deviation of the characteristic mass from the z = track 
is by more than a factor 2 at 2 ^ 20. At the highest metallic- 
ity considered in this analysis, [M/H] = —1, the character- 
istic mass in the line- (dust-) induced fragmentation is 100 
(0.1), 350 (0.7), 640 (3.3) M Q at z = 0,5,10, respectively. 
At 2 ^ 20, it is no longer possible to discriminate between 
the two fragmentation phases (this is why the long-dashed 
and dot-dashed lines in the upper branch are interrupted) 
and the fragmentation masses are ps IOMq. 

The above results are therefore consistent with the gen- 
eral trend suggested by Smith et al. (2009) , that the charac- 
teristic fragment mass Mf rag first is high, then it decreases, 
and then it is high again (3 modes) but in our analysis this 
trend is much less pronounced. To better compare with their 
analysis, we have run a set of test models with no dust and 
the characteristic masses in the line-induced fragmentation 
branch as a function of metallicity and redshift are shown 
in the right panel of Fig. [3] 

Since the total metallicity is the same in models with 
and without dust, about a factor of three more metals are 



present in the gas-phase in the no-dust cases. At z = 0, 
Mf ra g « 1000M Q when the metallicity is [M/H] < -4 
and Mfr a g w IOMq above this threshold. At z ^ 10, the 
CMB starts to affect the characteristic mass when [M/H] 
—3.5, increasing it to several 10s (100s) Mq when [M/H] 
= —3 (^ —2). At higher redshifts, the characterstic frag- 
ment masses are » 1000M Q when [M/H] < -4 and [M/H] 
> — 2 but never decrease below 100 Mq in the intermediate 
range —4 ^ [M/H] ^ —2. Therefore, our analysis confirm 
that the CMB inhibits fragmentation at large metallicities 
but, in the absence of dust grains, we do not find any metal- 
licity /redshift range where fragments of a few solar masses 
can form. As already pointed out in section 1, in their sim- 
ulation Smith et al. (2009) do not consider H2 formation 
heating nor H2 formation on dust grains which may explain 
why they tend to find smaller fragment masses in the in- 
termediate metallicity range. Note also that their model in- 
cludes more atomic species than just C and O, such as Fe, 
Si, and S, which provide additional cooling. 



4 CONCLUSIONS 

In this paper we have investigated the dependence of ther- 
mal and fragmentation properties of star-forming clouds on 
the environmental conditions, focussing on the interplay be- 
tween metals, dust grains and the CMB at different red- 
shifts. 

Our results, based on an extended grid of models with 
varying initial metallicity and redshift, can be summarized 
as follows: 

(i) In the absence of dust grains, when fragmentation oc- 
curs at the end of the line-cooling phase, moderate charac- 
teristic masses (of 10s of Mq) are formed at redshifts z ^ 10 
when the metallicity is 10~ 4 Zq ^ Z ^ 1Q~ 2 ' 5 Zq\ at lower 
metallicities, the clouds follow the primordial fragmentation 
mode and at higher metallicity, Z > 10~ 25 Z Q , the CMB 
inhibits fragmentation and only very large masses ~ 100s of 
Mq, are formed. These effects become even more dramatic 
at z > 10, when the fragmentation mass scales are always 
^ 100s of Mq, independent of the initial metallicity. 

(ii) When dust grains are present, at z = two regimes 
are present: the line-induced fragmentation mode, which 
leads to fragment masses > 50Mq at any metallicity, and the 
dust-induced fragmentation mode, which occurs at higher 
densities and leads to sub-solar mass fragments when Z ^ 
W~ e Z Q . The effects of the CMB at z > are impor- 
tant in both physical regimes: at very low metallicities, 
Z < 10~ 4 Zq, line-cooling becomes less efficient because of 
the suppression of H2 formation on dust grains when the 
grain temperature becomes ^ 75 K, at 2 > 20. At higher 
metallicities, metal line-cooling causes the gas to thermally 
couple with the CMB and the resulting fragments are very 
massive, with Mfrag > IOOMq for z > 5. Dust-cooling re- 
mains relatively insensitive to the presence of the CMB up to 
metallicities of Z ~ W~ 3 Zq and sub-solar mass fragments 
are formed at any redshift. Above this threshold, heating of 
dust grains by the CMB at 2 > 5 favors the formation of 
larger masses, which become super-solar when Z ^ W~ 2 Zq 
and 2 10. 

(iii) Clouds enriched by metals and dust grains with a 
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Figure 3. Characteristic fragmentation mass as a function of metallicity for 2 = (solid), 5 (dotted), 10 (short-dashed), 20 (long-dashed), 
and 30 (dot-dashed). The left panel shows the same set of models presented in Fig. [T] with the lower (upper) branch representing the 
masses associated to dust- (metal-) induced fragmentation. The right panel illustrates a test model with no dust. 



total metallicity Z ^ 1Q~ x Zq collapsing at z > 5, expe- 
rience only a single fragmentation epoch (line-induced and 
dust-induced phases merge) and the characteristic masses 
are controlled by the CMB temperature, with Mf rag = 
0.5, 3, 20M© for z = 10, 20, 30. 

The characteristic fragmentation masses are schemat- 
ically summarized in Fig. |4j where we report the typical 
scales associated to dust-induced fragmentation (left panel) 
and in the absence of dust (right panel) in the differ- 
ent redshift-metallicity regimes discussed above. We have 
not considered the fragmentation masses associated to line- 
cooling in the models with dust as they are always > 50Mq 
independent of the initial metallicity and redshift (see the 
upper branch in the left panel of Fig. [3J • 

Our results confirm the general trend suggested by pre- 
vious semi-analytic (Omukai et al. 2005) and numerical 
models (Smith et al. 2009), that the CMB tends to inhibit 
fragmentation of relatively metal-enriched clouds collapsing 
at high redshift. However, the interplay between metals, dust 
and the CMB generates a complex dependence of the frag- 
mentation mass on environmental conditions. 

It is important to stress that in order to explore the 
observational consequences of such an environmental de- 
pendence, it would be necessary to infer the variation of 
the functional form of the IMF and not only of the char- 
acteristic mass. This issue has been recently addressed by 
Clark et al. (2009) who predict that, in star forming re- 
gions with Z > W~ 5 Zq, competitive accretion can provide 
a natural explanation for the power-law form of the IMF, 
which is observed to be remarkably uniform in a wide vari- 
ety of environments. The competitive accretion of gas from 
a common reservoir by a collection of protostellar cores re- 
quires that the characteristic fragmentation mass be much 
smaller than the mass of the cooling gas available. As a con- 
sequence, Clark et al. (2009) suggest that there may be a 



link between the CMB and the slope of the IMF and that, 
in the absence of dust, competitive accretion will not oc- 
curr in regions with metallicities Z ^ 10~ 3 5 Zq collapsing 
at high redshift whereas it will occurr at all redshifts and 
metallicities if dust cooling can operate. 

As we have mentioned in section 1, one of the main 
motivation for the present study was to assess whether the 
CMB could be a viable way to modulate the stellar IMF 
with redshift as suggested by recent analyses of the ele- 
mental abundances of VMP and CEMP stars (Tumlinson 
2007; Komiya et al. 2007; Heger & Woosley 2008). In Fig. [5] 
we show the redshift dependence of the fragmentation mass 
scales found in the present analysis selecting only the pa- 
rameter space which is compatible with fragment masses in 
the range Mf rag ^ IQMq in the presence of dust grains. As a 
comparison, the thin dotted line shows the empirical depen- 
dence of the characteristic stellar mass on redshift proposed 
by Tumlinson (2007, see eq. 1). It is important to stress that 
the relation between fragmentation mass scale and charac- 
teristic mass shaping the stellar IMF is not straighforward 
and several processes, such as competitive accretion or dy- 
namical encounters, can act to make it not monothonic. 
However, the inspection of Fig. [5] suggests that for all but 
the highest metallicity, the redshift dependence predicted by 
our analysis is generally shallower than what inferred from 
the empirical IMF. Moreover, observations indicate that the 
fraction of CEMP stars decreases with metallicity being 
100% at [Fe/H] < -4.5, ~ 40% at [Fe/H] < -3.5 and ~ 20% 
at [Fe/H] < -2 (Beers & Christlieb 2005). If CEMP stars 
reflect the underlying IMF in the range 1M < M < 8M , 
our results seem to indicate an opposite trend: the effects of 
the CMB on the characteristic mass increase with metallic- 
ity, suggesting that the CMB can not be the only physical 
driver for the required modulation of the stellar character- 
istic mass. 
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Figure 4. Schematic represtentation of the characteristic masses associated to the relevant fragmentation regimes in the metallicity- 
redshift plane. In the left panel, we report only the typical scales associated to dust-induced fragmentation, since line-induced fragmenta- 
tion leads to masses > 50Mq independent of the initial metallicity or redshift. In the right panel, we show the scales expected in models 
with no dust (see text). 
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Figure 5. Redshift evolution of the fragmentation mass Mf rag 
for clouds with initial metallicities [M/H] = —4, —3, —2, —1 (from 
bottom to top). The thin dotted line shows the empirical depen- 
dence of the stellar characteristic mass on the CMB temperature 
proposed by Tumlinson (2007, see text). 
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